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O ■ ABSTRACT 

The prominent broad Fe II emission blends in the spectra of active galactic nuclei have been shown to vary 
in response to continuum variations, but past attempts to measure the reverberation lag time of the optical 
Fe II lines have met with only limited success. Here we report the detection of Fe II reverberation in two 
Seyfert 1 galaxies, NGC 4593 and Mrk 1511, based on data from a program carried out at Lick Observatory 
vQ ' in Spring 2011. Light curves for emission lines including Hf3 and Fe II were measured by applying a fitting 

routine to decompose the spectra into several continuum and emission-line components, and we use cross- 
correlation techniques to determine the reverberation lags of the emission lines relative to V-band light curves. 
In both cases the measured lag (rcen) of Fe II is longer than that of H/3, although the inferred lags are somewhat 
se nsitive to the choice of Fe II template used in the fit. For spectral decompositions done using the Fe II template 
of IVeron-Cettv et all (l2004i) . we find reen(Fe II)/rcen(H/3) = 1 .9 ± 0.6 in NGC 4593 and 1 .5 ± 0.3 in Mrk 15 1 1 . 
The detection of highly correlated variations between Fe II and continuum emission demonstrates that the Fe II 
,S^ ' emission in these galaxies originates in photoionized gas, located predominantly in the outer portion of the 

Qh, broad-line region. 

Q ' Subject headings: galaxies: active — galaxies: individual (NGC 4593, Mrk 1511) — galaxies: nuclei 
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1. INTRODUCTION 

Blends of Fe II emission lines are often among the most 
prominent broad emission features in the ultraviolet (UV) and 
optical spectra of broad-lined active galactic nuclei (AGNs), 
and the integrated flux of Fe II emission in quasars can be 
greate r than that of any other single emission line, including 
Ly a (IWills et al.lll985l) . Despite more than three decades of 
observational and theoretical effort, the physical conditions 
that give rise to Fe II emission have remained very difficult 
to determine. This is due in part to the complex energy-level 
structure of the Fe"^ ion and the very large number of indi- 
vidual Fe II emission lines that appear in broad blends. Re- 
cent spectral synthesis models for the Fe^ ion include hun- 
dreds of energy levels and tens of thousands of individual 
ti-ansitions (iVemeret al.lll999USigut & Pradhanll2003l) . A va- 
riety of processes in the broad-line region (BLR) can con- 
tribute to Fe II line production, including colli sional excitation 
as well as c ontinuum and hne fluorescence (iPhiUipsI 1 1978"bt 
iNetzer & W ills 1983| IComn-Souffrin et all 1198011 11 is not 
yet fully determined whether the Fe II production in AGNs 
occurs in gas heated solely by photoionization, or whether 
collisional ionization mi ght play a significant or even domi- 
nant role (I Colhn & Jolvir2000l: [Baldwin et ani2004l) . A par- 
ticularly promising recent development has been the realiza- 
tion that anisotropy in the emission from high column density 
clouds can have a strong impact on the observed Fe II spec- 
trum, and models for the emission from the shielded side of 
photoionized clouds have sh own significant impr ovements in 
fitting observed Fe II spectra dFerland et al.ll2009h . 

The Fe II emission lines are closely linked to several fun- 
damental issues in AGN physics and phenomenology. From 
principal-component analysis of quasar spectra, the relative 
strength of Fe II em ission is one of the major characteristics 
of "Eigenvector 1" (iBoroson & Greenl [l992). the component 
which accounts for much of the variance among the quasar 
population. Empirically, Eigenvector 1 represents the anti- 
coiTelations between Fe II and [O III] equivalent widths, and 
between Fe II strength and H/3 width. The fundamental physi- 
cal driv er of the E igenvector 1 sequence is th e Eddington ratio 
L/LEdd (ISulentic et al...20()0t lBorosonll2002h . but the physical 
mechanism responsible for the increase in Fe II strength with 
Eddington ratio is not fully understood. It may be linked to 
the presence of a soft X-ray excess at high Eddington ratio, 
which would produce a larger partially ionized zone of warm, 
Fe Il-emitting gas (e.g.. Marziani et al.i :2001: Boroson 2002), 
while o rientation cou l d play a secondary role (Marziani et al. 
1200 Ih . [Shields et al.l (120101) argue against the X-ray excess 
as the primary driver of Fe II strength, proposing instead that 
differences in gas-phase iron abundance, driven by selective 
depletion of iron onto grains, are piimaiily responsible for 
the wide range of observed Fe II li ne strengths in AGN s. In 
another interpretat ion proposed by iFerland et al.l (12009 ) and 
iDong et al] (1201 Ih . the competing forces of gravity and ra- 
diation pressure set a critical column density for clouds to 
remain gravitationally bound within the BLR, and at higher 
Eddington ratio, the higher column density of surviving BLR 
clouds would produce a larger Fe II/H/3 flux ratio. For AGNs 
having extremely strong Fe II emission, it has also been sug- 
gested that the shock-heated gas associated with circumnu- 
clear star formation m ay be responsible for the Fe II enhance- 
ment ( Lipari et al.llT993D . 

Furthermore, Fe II emission is potentially valuable as a 
tracer of chemical evolution. The strength of the UV Fe II 



emission relative to Mg II in quasar spectra has been used 
as a proxy for the iron to g-element abundance ratio in the 
BLR (e.g.,'Yoshii et al.l l 19981: [Dietrich et allllOOllKurk et all 
2007; Jiang et al. 200^, although the sensitivity of Fe II/Mg II 
line ratios to factors such as gas density and microturbulence 
means that observed line ratios are not st raightforward indi- 
cators of the u nderlying abundance ratios Cyemer et al.ll2003l : 
[Baldwin et all |2004 iBruhweiler & Verneij 12008b . Improved 
understanding of the physical conditions responsible for Fe II 
emission could have important ramifications for elucidating 
the history of metal enrichment in the densest regions of the 
Universe at high redshift. 

A key quantity of interest is the spatial scale of the Fe II- 
emitting region. However, there are very few measur e ments 
available that directly constrain its size. iMaoz et aTl (119931) 
carried out reverberation mapping of the UV Fe II fines in 
NGC 5548 and found a lag of about 10 days, similar to that 
of Lya, indicating that the UV Fe II lines originate from 
within the BLR. This has been the only measurement of a re- 
verberation lag for the UV Fe II lines to date. Aside from 
reverberation mapping, one additional direct constraint ex- 
ist s: a spectroscopic microlensing study of a lensed quasar 
bv iSlus e et al. (2007) found evidence that the Fe II emission 
originates largely in the outer portion of the BLR. 

Previous attempts at reverberation mapping of the optical 
Fe II blends have generally not led to clear detections of 
reverberation la gs. For the weU-studied AGN NGC 5548, 
iVesterg aard & P etersonI (|2b05) examined 13 years of moni- 
toring data and found that the Fe II flux responded to con- 
tinuum changes on timescales of less than several weeks, but 
the available data did no t allow for a definite measurement 
of the lag. iKuehn et alj ( l2008 f) measured light curves for 
the Fe II blends in Ark 120 and found a long-term variabil- 
ity trend that followed the continuum changes, but the data 
did not yield a significant cross-correlation lag. They con- 
cluded that the Fe II emission region in Ark 120 was likely 
to be several times larger than the H/3-emitting zone of the 
BLR. An alternative possibility was that the Fe II emission 
might be powered by collisional excitation rather than by pho- 
toionization, although they were unable to fit the Fe II emis- 
sion blends well with collisional excitation models. Bia n et alJ 
(2010) measured the F e II light curve for PG 1700H-518 using 
monitoring data from iKaspi et al] (120001) . While they were 
able to detect evidence of reverberation in the Fe II lines, 
the cross-correlation analysis did not yield a highly signif- 
icant peak, and the measured lag of 209^)47 days was very 
uncertain. Variability of the optical Fe II emission has been 
examined in other AG Ns (e.g., iG iannuzzo & Stirpe 199ot 
iDoroshenkoetal.111999^, iKollatschn y et al..i200Q : .Wang et all 
2005; Shapovalova et al. 2012), but most datasets have not 
been suitable for measurement of reverberation lags. The 
overall picture that emerges from these studies is that while 
optical Fe II emission does respond to continuum variations 
at least over long timescales, it tends to have a lower ampli- 
tude of variability than H/3, and it does not generally show a 
clear reverberation lag signature. 

In the absence of direct constraints on the relative sizes 
of the Fe II and H/3 emission regions, indirect clues have 
come from line-profile m easurements. Phillips (1978a ) and 
IBoroson & Green! (119921) found that the velocity widths of the 
Fe II lines are generally similar to that of H/3, suggesting that 
Fe II and H/3 originate from the same region within the BLR. 
Subsequent work has uncovered subtle systematic differences 
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between H/3 and optical Fe II widths, which hint at an origin 
for Fe II in the outer portion of the BLR or perhaps within 
an "intermediate-line region" corresponding to the transition 
between the BLR and the dusty toru s (Marziani et al. 2003; 
iPopovic et ai] |200?, 'Hu et al."2008b'; "Kovac evic et al.ll2010l) . 
For a sa mple of 40 00 Sloan Digital Sky Survey (SDSS) quasar 
spectra, iHu et alj (|2008b) found that the full- width at half 
maximum intensity (FWHM) of the optical Fe II lines is typ- 
ically about 0.75 X FWHM(H/3), albeit with substantial scat- 
ter, implying a size for the Fe Il-emitting region that is typi- 
cally about twice as large as the H/?-emitting zone of the BLR. 
Here, we present new Fe II reverberation-mapping results 
for two nearby Seyfert galaxies, NGC 4593 and Mrk 1511, 
from the Lick AGN Monitoring Project 2011. Section 2 gives 
a brief overview of our observing campaign. Section 3 de- 
scribes the fitting method used to isolate the continuum and 
emission-line components in the spectroscopic data and the 
measurement of light curves. In Section 4 we describe the 
cross-correlation lag measurements and the dependence of the 
lags on the choice of Fe II template, and Section 5 presents 
a discussion of the results. In an Appendix, we present a 
brief discussion of results based on additional newly-available 
Fe II templates. While the reverberation lag measurements 
are modestly sensitive to the choice of Fe II template, both of 
these objects clearly show evidence of Fe II reverberation in 
response to continuum variations. 

2. OBSERVATIONS AND REDUCTIONS 

Initial results from the Lick AGN Monitoring Project 
2011 were previously presented by Barth et al. (2011b) and 
iPancoast et alj (120121) . and full details of the spectroscopic 
and photometric observing campaigns, analysis procedures, 
and light-curve data for all targets will be presented in a forth- 
coming series of papers. Here, we briefly review the key as- 
pects of the observing program. 

2.1. Photometry 

From early March through mid- June 2011, we obtained 
queue-scheduled V-band images of NGC 4593 and Mrk 1511 
using the 0.76 m Katzman Automatic Imaging Telescope 
(KAIT) at Lick Observatory (Filippenko e t alj I2001h . flie 
0.9 m telescope at the Brigham Young West Mountain Obser- 
vatory (WMO), the Faulkes Telescope South at S iding Spring 
Observatory, and the Palomar 1.5 m telescope fcenko et al.l 
I2OO6). Our goal was to obtain nightly imaging for each tar- 
get, but weather and scheduling issues led to some gaps in 
coverage, particularly during the initial portion of the cam- 
paign. The temporal sampling cadence of our observations 
is described in Table [T] Exposure times were typically 180- 
300 s. All images were bias-subtracted and flattened, and 
cosmic-ray hits were removed using the LA-COSMIC routine 
(|yanDokkum2001). 

Image-subtraction photometry was carried out using a ver- 
sion of the ISIS code ( Alard & Lupton.1998.) rnodifie d by 
the High-z Supernova Search Team (ITonrv et al.l 120031) and 
with additional modifications by W. Li, and also using the 
HOTPANTS package by A. BeckeQ which is based on the 
methods described by lAlardI (120001) . For each telescope, a 
high-quality template image was chosen, and the template 
was then aligned with each night's image and convolved with 
a spatially varying kernel to match the point-spread function 
of that image. After subtracting the scaled template image, the 

^"^ http://www.astro.washington.edu/users/becker/c_software.html 



variable AGN flux is left as a point source in the subtracted 
image, allowing for aperture photometry using the IRAF0 
DAOPHOT package. The photometric aperture radius used 
for each telescope was set to match the average point-source 
FWHM for images from that telescope. 

Light curves were initially constructed separately with the 
data from each telescope, and then combined together to as- 
semble a single light curve for each AGN. Normalization was 
done by modeling the variability using Gaussian processes 
to obtain a finely-sampled model version of the light curve 
for each AGN , following the same procedure described by 
IPancoast et alj (jJOlli). Additive and multiplicative scaling 
factors were applied to match the modeled light curve to each 
telescope light curve, using the WMO light curve as the ref- 
erence data for the model. The normalization code simulta- 
neously constrained all scaling factors to obtain the best fit 
between all the telescope light curves using a Markov Chain 
Monte Carlo algorithm. Then, scaling factors determined 
from the modeling were applied to the image-subtraction light 
curves, to normalize the data from all telescopes to a com- 
mon flux scale. The final light curve was prepared by taking a 
weighted average of observations taken within 6 hours of one 
another. The V-band light curves are shown in Figure [T] in 
arbitrary units of differential flux. 

2.2. Spectroscopy 

Our Spring 201 1 observing campaign consisted of 69 nights 
at the Lick 3-m telescope between 201 1 March 27 and June 13 
UT. Observations were done using the Kast dual spectrograph 
dMiller & Stonel[T993b . A D55 dichroic was used to separate 
the blue and red beams with a crossover wavelength of about 
5500 A. In this paper, we discuss only measurements from 
the blue arm of the spectrograph, where we used a 600 lines 
mm"' grism over ^ 3440-5520 A at a scale of 1.0 A pixel"'. 
All observations of these two AGNs were done with a 4"-wide 
slit oriented at a position angle of 45°. Standard calibration 
frames including arc lamps and dome flats were obtained each 
afternoon, and flux standards were observed during twilight. 
Exposure times were normally 2 x 600 s and 2 x 900 s for 
NGC 4593 and Mrk 1511, respectively. The weather at Mt. 
Hamilton was somewhat worse than average during Spring 
2011; we were able to observe NGC 4593 on 43 of the 69 
nights, and Mrk 1511 on 40 nights. 

Spectroscopic reductions and calibrations followed stan- 
dard methods implemented in IRAF and IDL. A large ex- 
traction width of 10"3 was used in order to accommodate 
the full extent of the AGN spatial profiles observed on nights 
with poor seeing. Error spectra were extracted and propa- 
gated through the full sequence of calibrations. In the reduced 
spectra, the median signal-to-noise ratio (S/N) per pixel in the 
range 4600-4700 A is 122 for NGC 4593 and 71 for Mrk 
1511. 

3. SPECTROSCOPIC DATA ANALYSIS 

3.1. Spectral Fitting Method 

The reduced spectra were first normalized to a 
uniform flux scale by emplo ying the procedure of 
Ivan Groningen & WandersI (119921) . This method applies 

^' IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in Astron- 
omy, Inc., under cooperative agreement with the National Science Foundation 
(NSF). 
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Fig. 1. — Light curves of the V band, the AGN featureless continuum (FC), H/3, H7, He II (for Mrk 1511 only), and Fe II. The Balmer-line and He II light 
curves contain flux from narrow components, including [O III] A4363 in the H7 light curve. The V-band photometric light curves are given in differential flux 
units (i.e., relative to the flux in a reference image), as determined by t he im age-subtraction photometry procedure. The FC and emission-line light curves are 
those measured from the spectral fits which used the lVeron-Cettv et alj 120041) Fe II template. 



a flux scaling factor, a linear wavelength shift, and a Gaussian 
convolution to each spectrum in order to minimize the 
residuals between the data and a reference spectrum con- 
structed from several of the best-quality nights. The scaling 
is determined using a wavelength range containing the [O III] 
A5007 line, which is assumed to have constant flux. The 
upper panels of Figure |2] show the mean of all scaled spectra 
for each object. 

To assess the accuracy of the scaling procedure, we mea- 
sured the light curve of the [O III] A5007 line in the scaled 
spectra and calculated the residual "excess" scatter a^ in the 
[O III] light curve. We fo llow the usual defin ition for the nor- 
malized excess variance JNandra et al.lll99'7b : 



^'-w^^^^^-- 



^^r 



(1) 



where A^ is the total number of observations, /i is the mean 
flux, and X, and cr, are the individual flux values and their un- 
certainties. Then cr^ gives a measure of the fractional scat- 
ter in the light curve due to random errors in flux scaling, 
over and above the scatter that would be expected based on 
the propagated photon-counting uncertainties. For NGC 4593 



and Mrk 1511, we obtain an excess scatter of cr^ = 1.2% and 
2.0%, respectively. These low values indicate that the scaling 
procedure worked well. As we describe below, these mea- 
sures of residual scatter will be added to the error budget of 
the emission-line light curves to give a more realistic estimate 
of the total uncertainties. 

Traditionally, broad emission-line fluxes in reverberation- 
mapping data have been measured by choosing continuum 
regions on either side of an emission line, fitting a line to 
the continu um, and integrating the flux above that line (e.g., 
iKaspi et 311 :2000). However, in AGN spectra there may be 
no pure continuum regions at all surrounding H/?, because 
Fe II, He II, and other emission features are present. Choosing 
continuum regions that contain some emission-line flux could 
potentially bias a reverberation lag measurement, particularly 
for velocity-resolved reverberation signals in the faint high- 
velocity wings of H/3. An alternative approach is based on 
decomposition of spectra into continuum and emission-line 
components, so that the flux contributions of individual emis- 
sion features can be isolated. Our mult i-component fitting 
method was previously applied to Mrk 50 ('Barth et al."2 01 Ibl) 
and is similar to the method used by Park et al. (2012). We 
carry out fits over the wavelength range extending from just 
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TABLE 1 
Light-Curve Sampling Statistics 



Dataset 


A'obs 


AWdian (days) 


Armean (days) 


NGC 4593 Photometry 


75 


0.95 


1.12 


NGC 4593 Photometry (Averaged) 


57 


1.01 


1.48 


NGC 4593 Spectroscopy 


43 


1.02 


1.81 


Mrk 1511 Photometry 


71 


1.00 


1.51 


Mrk 1511 Photometry (Averaged) 


59 


1.05 


1.82 


Mrk 1511 Spectroscopy 


40 


1.01 


1.94 



Note. — Wq^s is the number of observations in each light curve, and A^njeji;,,, 
and A?mcaii give the sampling cadences, defined as the median and mean separa- 
tion between adjacent data points in each light curve. Photometric light curves 
were condensed by taking a weighted average of any observations separated by 
less than 6 hours. 
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Fig. 2. — Mean and rms spectra of Mrk 1511 andNGC4593. Upper panels: Mean spectrum (black) overlaid with the model fit done using the lVeron-Cettv et aP 
1200 41) Fe II template (red). Individual fit components include the AGN power law (green), old stellar population (puiple), Fe II emission (grey), H/3 (magenta), 
[0 III] (blue). He II (fight blue), and He I (orange). Lower panels: The standard iTns spectrum (top) and an rms spectram constructed from individual nightly 
spectra after the AGN power-law and stellar components were removed (bottom). 



bluewai-d of the H7 line (Are.st « 4150 A) up to 5470 A near 
the dichroic cutoff. Over this wavelength range, each nightly 
spectrum was fit with a model consisting of several compo- 
nents: a power-law featureless continuum (PC), a starlight 
template broadened in velocity by convolution with a Gaus- 
sian, emission Unes including H;3, [O III] AA4959,5007, He II 
A4686, and three He I Hnes (4471, 4922, and 5016 A), and 
an Fe II template broadened in velocity by convolution with a 
Gaussian. 

For the starlight template, we use an 11 Gyr-old sim- 
ple stellar population m odel at solar metallicity from 
iBruzual & Charlofl (J2003). We experimented with adding 
younger stellar population components, but the 1 1 Gyr model 
provided a sufficiently good fit that more complex stellar pop- 
ulation models are not warranted, and the flux of a young pop- 
ulation component is usually not very well constrained over 
this limited wavelength range. 

We used two different templates to model the Fe II lines, 
from Boroson & Green (1992; hereinafter BG92), and Veron- 
Cetty et al. (2004; hereinafter V04). The overall quaUty of 
the fits is very similar with the two Fe II templates, but there 



are some subtle differences that affect the reverberation lag 
measurements. We find that the V04 template produces bet- 
ter results than the BG92 template (see Section 14.4) , but we 
describe measurements and reverberation lags derived from 
both sets of spectral fits in order to illustrate the dependence 
of the results on the choice of Fe II template. 

The broad component of H/3, as well as the [O III] lines, are 
represented by fourth-order Gauss-Hermite functions, while 
Gaussians were used to represent weaker features including 
narrow H/3, the broad and narrow components of He II, and 
the He I lines. We found that the complex shape of the 
H7H-[0 III] A4363 blend was difficult to model, and instead of 
adding more model components we simply set the wavelength 
range subtended by this feature to have zero weight in the 
fit (4280^395 A rest wavelength). Additionally, a uniform 
foregro und extinction was applied to the model spectrum us- 
ing the ICardelli et al.l ([1989) reddening law. This extinction 
correction represents the combined effects of Galactic fore- 
ground extinction, extinction within the AGN host galaxy, and 
wavelength-depende nt slit losses resul ting from the nonparal- 
lactic slit orientation (iFilippenkol 1 982h . Since these slit losses 
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Fig. 3. — Same as Figure|2] but for fits done using the lBoroson & GreenI j 19921) Fe II template. 

differ from night to night, we allowed the value of E(B-V) to 
vary freely in each fit. 

The full model includes 29 free parameters, and was fit- 
ted to the observed spectrum by y- minimization using a 
Levenberg-Marquardt technique (Markwardt 2009). For each 
galaxy, the fit was first carried out on the high-S/N mean spec- 
trum, and then the parameters from the best fit to the mean 
spectrum were used as starting parameter estimates for the fit 
to each nightly spectrum. 

For both objects, the fit optimization always drove the flux 
of He I A4471 to zero, using either Fe II template. However, 
inclusion of the other He I Unes at 4922 and 5016 A does 
significantly improve the fit quality compared to model fits 
that omit these lines, at least when using the V04 template. 
Previous work has shown that these two He I lines can be im- 
portant contributor s to the "red-shelf" region redward of H/3 
dVeron et al.ll2002h. Th e red shelf also contains some Fe II 
emission ( lKoristalll992l) and the I Zw 1 templates contain 
some Fe II flux in this region. However, the relative amount 
differs between the BG92 and V04 templates, such that the 
V04 template has relatively weaker Fe II emission in the H/3 
red shelf, and this leads to noticeably different fitting results 
for the two templates. Using the V04 template, we found that 
the Fe II template alone did not produce an adequate fit to the 
red shelf unless the He I lines were added as separate compo- 
nents. When fits are carried out using the BG92 Fe II template 
instead, the fluxes of the He I A A4922, 5016 components in the 
H/3 red shelf region go to zero or nearly zero, and essentially 
all of the flux in the H/3 red shelf is taken up by the Fe II 
template. 

Another difficulty in fitting this region is that the He I A4922 
line is degenerate with the red wing of H/3; in order to achieve 
a consistent deblending of this spectral region into H/3 and 
He I flux, we constrained t he 4922 A and 5016 A lines to have 
identical fluxes (following lVestergaard & Petersonll2005h and 
identical velocity widths. In the fits, these two components 
essentially track the same shape as two broad bumps in the 
broadened Fe II templates. Due to the degeneracy between the 
several features contributing to the H/3 red shelf, we are un- 



able to determine a physically unique decomposition of this 
region into separate contributions from H/3, Fe II, He I, and 
[O III]. The fits are able to match the observed spectra well 
in this region using either Fe II template, but the slight differ- 
ences in the red shelf decompositions may be partly responsi- 
ble for the dependence of the reverberation lags on the choice 
of template. 

The strength and broadening of the Fe II component in the 
fit are primarily determined by the prominent blend around 
Aiest = 4500-4700 A containing several transitions from Fe II 
multiplets 37 and 38. The velocity-broadening width of Fe II 
was allowed to vary independently to optimize the fit to each 
nightly spectrum, but the variation from night to night was 
relatively small. From the best fit to each night's data us- 
ing the V04 template, the median velocity broadening (i.e., 
the dispersion of the Gaussian broadening kernel) applied to 
the Fe II template is 1885 ±75 km s"' for NGC 4593 and 
1635 ±50 km s"' for Mrk 1511. The V04 template itself 
was constructed using Lorentzian profiles for each line, with 
a width of FWHM =1100 km s"' to match the line profiles 
in I Zw 1 ( Veron-Cetty et al. 2004) . We find a similarly small 
night-to-night scatter in the width of the Fe II broadening ker- 
nel when using the BG92 template (±65 km s"'), although 
the overall amount of broadening is substantially smaller due 
to the different template structure. 

The BG92 template is characterized by an intrinsic 
FWHM«900 km s"' (e.g., iHu et alj l2008bl) . but as an 
empirically-derived template it contains features with a range 
of widths, and shapes that do not correspond precisely to sim- 
ple Gaussian or Lorentzian profiles. As a result, it is not 
straightforward to compare the inferred Fe II widths as mea- 
sured using the two different templates, but for each template 
the fitting procedure consistently converges on a value for the 
Fe II broadening that varies only by a few percent from night 
to night. We comment further on the Fe II profile widths in 
Section|5] 

Figures |2] and |3] show the mean spectrum for each AGN 
with its best-fitting model, and the individual fit compo- 
nents, for the fits done using the V04 and BG92 Fe II tem- 
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plates. The root-mean-square (rms) spectrum is constructed 
by taking the sta ndard deviation of the scaled spectra (e.g., 
iKaspi et aI]|200Qi) . In each of these figures we show the rms 
spectrum, as well as a modified version constructed by sub- 
tracting the AGN EC and stellar-continuum components from 
each nightly spectrum before calculating the rms. Removing 
the continuum components gives an improved rms spectrum 
by eliminating sm all residual stellar absorption-line features 
JParketal .1120 121) . and the continuum-subtracted rms spec- 
trum more accurately depicts the rms variability profiles of the 
emission lines. Fits done using the BG92 template assign rel- 
atively more flux to Fe II than fits with the V04 template. For 
the BG92 fits, the increased Fe II flux (at the expense of the 
EC component) results in a higher rms flux in the continuum- 
subtracted rms spectra. 

3.2. Light-Curve Measurements 

Light curves were measured for individual emission 
lines from continuum-subtracted spectra. For H/3 and the 
H7H-[0 III] blend, we produced a residual spectrum by sub- 
tracting all other model components, and then measured the 
emission-line light curves by summation of the residual flux. 
The H7 and H/3 integration regions were 4320^450 A and 
4800-5000 A for NGC 4593, and 4420^560 A and 4920- 
5120 A for Mrk 1511 (in the observed frame). The resulting 
light curves include the constant contributions of narrow H/3 
and H7 as well as [O III] A4363 blended with H7. 

Light curves for He II A4686 were measured by summation 
of the best-fitting broad and narrow He II model components. 
This produced a less noisy light curve than the alternate ap- 
proach of calculating a He II residual spectrum for each night 
by subtracting all of the other model components. However, 
for NGC 4593 the He II emission was too weak to produce a 
useful light curve. The Fe II light curves were computed by 
summation of the best-fitting Fe II model over 4400^900 A, 
but any wavelength range would yield the same overall light- 
curve shape since the model fit assumes a uniform flux scaling 
factor for the Fe II template. The EC light curves were mea- 
sured by integration over the same wavelength range used for 
Fell. 

As described above, the residual scatter measured from the 
[O III] light curves represents a major contribution to the error 
budget in the light-curve measurements. Therefore, we com- 
bine the residual scatter (1.2% for NGC 4593 and 2.0% for 
Mrk 151 1) in quadrature with the measurement uncertainties 
in the individual light-curve points to obtain a more realistic 
estimate of the true overall uncertainties. The emission-line 
light curves displayed in Figure [T] are shown with these ex- 
panded error bars. Addition of the residual flux-scaling scat- 
ter to the error budget widens the measurement uncertainties 
in the cross-correlation lags described below, but only slightly. 
In Figure [U we show only the spectroscopic light curves mea- 
sured using the V04 template fits because the light curves 
measured from the BG92 fits are extremely similar in over- 
all appearance. 

Observations of other AGNs have shown that the amplitude 
of Fe II variability is typically somewhat smaller than that 
of H/3. For example, in NGC 5548 Vestergaard & Petersonl 
(12005) found that the Fe II variation amplitude was '^ 50%- 
75% of that of H/3. To characterize the variability amplitude 
in NGC 4593 and Mrk 1511, we use a-^, as defined in Equa- 
tion 1, which gives a measure of the normalized rms variabil- 
ity, corrected for measurement uncertainties. Based on the 



V04 fits, the results for NGC 4593 are (jAW)= 0.199 and 
crx(Fe II)= 0.138, while for Mrk 1511 we measure crx(H/3)= 
0.102 and crx(Fe II)= 0.085. With the BG92 template fits, the 
a^ values for H/3 are essentially unchanged, while o'x(Fe II) 
increases by ^ 10% relative to the values determined from 
the V04 fits. Thus, both objects show a somewhat higher am- 
plitude of variabilit y in H/3 than Fe II, similar to t he behavior 
seen in NGC 5548 /Vestergaard & Petersonll2005l) . 

For AGNs having a very low starlight fraction, it is of- 
ten possible to measure the AGN continuum flux from the 
spectra directly, without carrying out a decomposition to sep- 
arate the AGN fl ux from the host galaxy contribution (e.g., 
iKaspi et ani2000) . However, this approach is much less suc- 
cessful for AGNs having a substantial host galaxy component 
in the spectra, since the starlight dilutes the variability of the 
AGN continuum and introduces additional random noise due 
to nightl y variations i n seeing and target centering within the 
slit (e.g., iBentz et al.' 2008). Given the substantial contribu- 
tions of both stellar continuum and Fe II in these two objects, 
the spectral decomposition method is necessary in order to 
measure the AGN continuum flux without contamination or 
dilution by these other components. 

4. MEASUREMENT OF REVERBERATION LAGS 
4.1. Methods 

In order to measure the cross-correlation function (CCF) 
for unevenly sampled time series, we employ the in- 
terpolation cross-correlation function (ICCF) methodology 
and Monte Carlo e rror-a n alysis techniques describe d by 
Gaskell & Peterson ( 19871) . IWhite & Petersonl (119941) . and 
iPeterson et al.i (i2004i) . Emission-line lags for H/3, H7, He II, 
and Fe II were measured relative to both the V-band and EC 
light curves. CCEs were computed over a temporal range of 
-20 to +40 days in increments of 0.25 days. For each CCF we 
compute two measures of the lag time: Tpeak, which is the lag 
at the peak of the CCF, and Teen, the ce ntroid of the C CF for 
all points above 80% of the peak value (IPeterson et al.l l2004). 
Table|2]lists the lag values. The quantity r^ax listed in Table|2] 
gives the peak amplitude of the CCF and is a measure of the 
significance of the correlation between the two light curves. 
We obtain higher rmax values for NGC 4593 than for Mrk 
1511, primai-ily owing to the higher S/N of the NGC 4593 
light curves. Lag values are given in the observed frame and 
can be converted to the AGN rest frame by dividing by 1 +z, 
where the redshifts are z = 0.009 for NGC 4593 and z = 0.0339 
for Mrk 1511. 

There are advantages and drawbacks to using either the 
V-band data or spectroscopic EC light curves in the cross- 
correlation analysis. The V-band light curves have more fre- 
quent sampling than the spectroscopic data, and better tem- 
poral sampling will always improve the determination of the 
CCF. Additionally, the spectroscopic data suffer from system- 
atic effects including variable slit losses due to miscentering 
and differential atmospheric refraction, while the photomet- 
ric data are not susceptible to these problems. However, the 
V-band data include contributions of emission-line flux in ad- 
dition to the AGN continuum, while the spectroscopic decom- 
positions can be used to produce a "pure" EC light curve with- 
out contamination by strong emission lines. A disadvantage 
of using the spectroscopic data in cross-correlation measure- 
ments is the effect of correlated flux-calibration errors be- 
tween the emission-line and continuum light curves, which 
can introduce a spurious signal at zero lag. We find that using 
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Fig. 4. — Cross-correlation functions for H/3, H7, He II, and Fe II relative 
to the V-band photometric light curves. The autocorrelation function of the 
V-band light curve is shown in grey. 

the V-band light curves produces consistently higher-quality 
results, but the reverberation lags are consistent when mea- 
sured against either the V band or FC light curves. For com- 
pleteness. Table |2]presents all of the lag measurements based 
on the spectral decompositions using both the V04 and BG92 
Fe II templates, and for cross-correlations computed with both 
the V-band and FC light curves. 

We tested the effect of detrending the data by subtracting a 
linear fit to the light curves before computing the CCFs. As 
described by Welsh tll999h . removing long-term secular vari- 
ations from light curves can often improve the accuracy and 
significance of cross-correlation results. For Mrk 1511, de- 
trending had almost no effect on the CCFs, because the light 
curves are not characterized by any long-term brightening or 
fading. NGC 4593 does show a net increase in luminosity 
over the duration of our monitoring program, but CCFs com- 
puted using detrended data were actually of poorer quality 
than the original CCFs, having values of rmax that were lower 
by ~ 0.1. Therefore, we chose not to use detrended light 
curves for our final measurements. 

4.2. Reverberation lag results 

We first describe the results based on light curves measured 
from the V04 template fits. Co mparison with BG92 template 
fits is presented in Section l4!4l Figure |4] illustrates the CCFs 
measured relative to the V-band data for emission-line light 
curves measured from the V04 template fits, as well as the 
autocorrelation function (ACF) of the V-band light curve. For 
both AGNs the Fe II CCF has a shape roughly similar to that 
of the Hf3 CCF, but systematically shifted toward longer lags. 

Normally, the ICCF method computes two versions of 
the CCF: first by interpolating the continuum light curve 
(shifted by each trial value of the lag) to the temporal steps 
of the emission-line light curve, and then by interpolating the 
emission-line light curve to the time sampling of the contin- 
uum light curve. With high-quality data, the two resulting 
CCFs are normally very nearly id entical, and the final CCF i s 
taken to be the mean of the two dGaskell & Petersonlll987l) . 
However, when the two light curves have significantly dif- 
ferent temporal sampling frequency, or when the light curves 
have poor S/N, the two interpolated CCFs can sometimes dif- 
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Fig. 5. — Cross-correlation functions for Fe II relative to the AGN feature- 
less continuum (FC) measured from the spectral fits done using the BG92 and 
V04 Fe II templates. 

fer significantly in shape, and occasionally one of the two 
CCFs may contain features that are clearly spurious while 
the other has a more regular appearance. In such a situa- 
tion it may be preferable to use just one of the two interpo- 
lated CCFs rather than taking their mean. For NGC 4593, the 
two interpolated CCFs are nearly identical, and we base our 
measurements on the mean CCFs as illustrated in Figure |4] 
However, for Mrk 1511, we found that the second version of 
the CCF (i.e., from interpolating the Fe II light curve to the 
time steps of the V band) had a large and unphysical "notch" 
appearing just at the CCF peak. The other interpolation pro- 
duced a smooth and regular CCF. Overall the two interpolated 
CCFs had nearly identical shapes and widths, and they show 
the same overall lag response aside from the notch appearing 
at the peak of the first CCF. The presence of this notch in the 
mean CCF rendered the measurement of Tpeak and Teen sus- 
pect. Since the first interpolated CCF results in a smoothly 
shaped peak, for Mrk 1511 we use only the first interpolated 
version of the CCF to calculate Tpeak and Teen, rather than tak- 
ing the mea n of the two interpolations. This CCF is shown in 
Figure 130 

The H/3 lags measured relative to the V-band light curves 
are in the range of ^ 3-4 days for NGC 4593 and ^ 5-7 days 
for Mrk 1511. Due to asymmetry in the CCFs, the values of 
Teen and Tpeak differ slightly, but not by a significant amount 
considering the uncertainties in the lag measurements. The 
observed trend of H/3 having a longer lag than H7 is expecte d 
based on previous measurements (e.g.. iBentz et al.l l2010al) . 
Similarly, the fact that the He II lag is unresolved in Mrk 1511 
is consistent with expectations based on previous measure- 
ments. In AGNs having H/3 lags of < 10 days, the He II lag 
time is usually undetectable in datasets w ith nightly sampHng 
(e.g.. lBentz et alJl2010atlBarth et al.ll20r ib). The Fe II emis- 
sion has a longer lag than H/3 for both AGNs. In both cases the 
cross-correlation between the Fe II and V-band light curves 

^* For Mrk 1511, if we use the mean of the two interpolated fornis of the 
CCF instead of the first interpolated version, we obtain Tccn = 7.67^2 '2 days. 
This is consistent with the value listed in Table|2]but with a larger uncertainty 
due to the peculiar shape of the peak region of the mean CCF. The value of 
Tccii(H/3) is nearly unchanged if we use the mean CCF instead of the first 
interpolated version. 
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Cross-Correlation Lag Measurements 




Measurement 
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Tccn (days) 
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V04 Template Fits 
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H7 vs. V 
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-, 59+0:95 
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H/3 vs. FC 
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Fe II vs. FC 


0.88 
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57S-i 
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Mrkl511: 
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-0-72^il 
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H7 vs. V 
H7 vs. FC 
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•3 47+0.70 
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■X O7+0.98 
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H/3 vs. V 
H/3 vs. FC 
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7 00+1-" 
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Fe U vs. V 


0.80 


8.63+-3 
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Fe II vs. FC 
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Note. — All lags are given in the observed frame. Lag values measured from cross- 
correlations between two spectroscopic light curves (i.e. , between emission lines and FC, 
or between two emission lines) are unreliable as a result of correlated errors, which intro- 
duce a spurious signal at zero lag. This bias particularly affects the Tpeaj^ values. Cross- 
correlations measured relative to the V-band Hght curve are not susceptible to this bias. As 
described in the text, we conclude that the cross-correlations done relative to the V-band 
light curves, done using the V04 template fits, are the most reliable. 



produces a significant CCF peak, ■with r^ax = 0.88 for NGC 
4593 and 0.80 for Mrk 1511. 

We also measured the emission-line lags relative to the 
AGN FC component light curves. The resulting Teen and 
Tpeak valucs are generally consistent, ■within the uncertainties, 
■with the lag values measured against the V-band data. How- 
ever, the CCFs measured using the FC light curves have lower 
peak amplitude rmax than the corresponding CCFs measured 
against the V band, which we attribute to the poorer temporal 
sampling and greater calibration uncertainties of the spectro- 
scopic data. Thus, we consider measurements relative to the V 
band to be the best determinations of the emission-line lags, 
while the lags measured relative to the FC serve as a useful 
consistency check. As another check, we measured the FC 
light curves over a different wavelength range (5050-5150 A 
rest wavelength) and measured the cross-correlations of the 
H/3 light curves relative to this revised FC light curve. The 
measured lags only changed by a negligible amount (much 
smaller than the 1 a uncertainties) in comparison with the lags 
measured to the FC light curves as listed in Table |2] 

Figure |5] shows the Fe II vs. FC CCFs for both sets of tem- 
plate fits. For Mrk 1511, the choice of Fe II template has very 
little impact on the CCF shape or the lag in this measurement. 
However, the difference between the two templates is signif- 
icant for NGC 4593, where the light curves measured from 
the BG92 template produce a very jagged and irregular CCF 
and yield very unreliable lag measurements with large uncer- 
tainties on Teen and Tpeak- The V04 template, on the other 
hand, gives a better-behaved result and smaller lag uncertain- 
ties. The cause of this difference is not obvious, but it does 
appear that the V04 template gives a clearly superior result at 
least for NGC 4593. 

An additional possible concern is degeneracy between the 
stellar and AGN continuum components in the fit. For both 
objects, the light curves of the starlight component do not 
show any significant time-dependent features other than ran- 



dom noise. We measured cross-correlations between the stel- 
lar continuum component and the V-band light curves in order 
to test whether the stellar continuum light curves might con- 
tain some residual AGN flux. The resulting CCFs had very 
low peak ampHtudes of r^ax = 0.31 and 0.43 for NGC 4593 
and Mrk 1511, respectively, and had such irregular structure 
that no clear lag signal could be measured in either case. This 
provides further evidence that the spectral decompositions are 
properly separating the contributions of the continuum com- 
ponents. 

4.3. Estimating biases in reverberation lags 

A small bias in reverberation lag measurements can occur 
when broad-band photometry is used for the continuum flux 
measurement, due to the contribution of emission-line flux to 
the photometric light curves. The passband of the V filter 
includes flux from both H/3 and Fe II, both of which lag the 
continuum variations. When cross-correlating spectroscopic 
emission-line light curves against the V-band light curves, this 
emission-line contamination in the V-band data will tend to 
bias the cross-correlation lagtoward values that are lower than 
the true lag (Earth et al. 20 1 lab . To determine the level of this 
contamination for Mrk 1511 as an example, we extrapolated 
our mean-spectrum fit over a broader wavelength range (up 
to 6500 A) and carried out synthetic V-band photometry on 
the individual model components. In the mean spectrum, the 
ratios of emission-line to AGN featureless continuum (FC) 
flux as measured through a V-band filter are /(H/3)//(FC) = 
0.073, and /(Fe II)//(FC) = 0.155. 

We simulated the effect of the emission-line flux in the 
V-band light curves by creating mock reverberation datasets 
having properties sim ilar to those of our obser ved data. Fol- 
lowing the method of iTimmer & Koni^ (Il995h . we began by 
generating simulated AGN continuum light curves on a finely 
sampled temporal grid with bin size equal to 0.01 day, based 
on a power-density spectrum of the form P{f) (x /"^'^ (sim- 
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ilar to slopes measure d from Kepler AGN light curves by 
iMushotzkv et al.ll201 Ih . Each light curve was generated over 
a total duration of 1000 days. A shorter segment of length 
100 days was then selected randomly and normalized to have 
an rms variability equal to 15% of the median flux. The H/3 
and Fe II lines were assumed to respond linearly to continuum 
variations, with a delta function as the transfer function for 
simplicity, and with lag times of 6 days for H/3 and 9 days for 
Fe II. Variability amplitudes for the emission lines were nor- 
malized so that H/3 had the same 15% rms variability as the 
continuum, while Fe II was set to have a lower rms variabil- 
ity amplitude of 12%. Then, a "contaminated" V-band light 
curve was created by adding scaled versions of the H/3 and 
Fe II light curves to the EC light curve, in the relative propor- 
tions listed above for Mrk 1511. All of the light curves (EC, 
V band, H/3, and Fe II) were then degraded to approximate 
the sampling and S/N of real data. From the finely sampled 
light curves, one flux point was chosen per night, randomly 
sampled from a window of ±2 hours relative to a uniform 
24-hour cadence, and over a total spectroscopic monitoring 
duration of 80 days. Weather losses were approximated by 
randomly removing 10% of the points from the EC and V- 
band light curves and 40% of the points from the emission- 
line light curves. Random Gaussian noise was added to attain 
S/N = 100 for the EC and V-band data, and S/N = 50 for the 
emission-line light curves. The Fe II light curve was cross- 
correlated against both the pure EC light curve and the con- 
taminated V-band light curve, to determine the typical level of 
the bias, and this procedure was repeated 10"* times to build 
up a distribution of lag measurements for different initial re- 
alizations of the light-curve shape. 

We compiled the values of Teen for the subset of simula- 
tions that produced "successful" measurements of lag, mean- 
ing that the CCF had r,nax > 0.6 and yielded a nonnegative 
value of Teen; by these criteria 83% of the simulations were 
successful. For this subset, the median lag of Fe II relative 
to the pure AGN continuum is Teen = 9.1 ± 1.8 days, which 
is closely consistent with the input lag of 9 days. (The un- 
certainty of ±1.8 days represents the 68% confidence inter- 
val on Teen for the set of successful simulations.) When the 
simulated Fe II light curves are cross-correlated against the 
simulated V-band data including emission-line contributions, 
we find a median lag of Teen = 8.4 ±2.1 days. As expected, 
this is shorter than the lag measured with respect to the pure 
AGN continuum, but the difference of 0.7 days in the median 
measurements is smaller than the uncertainty in the measured 
Fe II lag for either of our targets. Thus, the impact of this bias 
on the emission-line lag measurements is fairly small, most 
likely below the level of the la uncertainties in the Fe II lags. 
The effect of emission-line contamination in the V-band light 
curves would be similar for all of the cross-correlation mea- 
surements for a given AGN, so the inferred relative sizes of 
the Fe II and H/3 emission regions would be only modestly 
affected. 

4.4. Differences due to Fe II Template Structure 

Among the reverberation measurements listed in Table |2] 
the strongest differences between results based on the V04 
and BG92 template fits are in the lag of Fe II relative to the V 
band. The BG92 fits give Teen values that are 2-3 days lower 
than the values measured from the V04 fits. This difference 
appears to stem from the different proportions of FC and Fe II 
fluxes for the two sets of spectral fits. As can be seen in Fig- 
ures |2] and |3] the broadened V04 template approaches zero 



flux density for the Fe II emission in the blue wing of H/3, 
as well as at the blue end of our fitting region. The broad- 
ened BG92 template, on the other hand, does not reach zero 
flux density anywhere in our fitting region, producing a higher 
"pseudo-continuum" level from blended Fe II lines and forc- 
ing the FC component to a lower overall flux than in the V04 
fits. We suspect that the BG92 template fits are including a 
small amount of AGN continuum flux into the Fe II fit com- 
ponent. This small contribution of AGN continuum flux in 
the Fe II light curves would add a spurious signal at zero lag 
to the cross-correlations between Fe II and the V-band data, 
slightly biasing the lags toward low values. 

These differences stem from the methods used to construct 
the Fe II templates. The BG92 template is based on an ob- 
served spectrum of 1 Zw 1 . Emission lines from several tran- 
sitions other than Fe II were removed, and the continuum was 
removed by fitting a polynomial to regions between strong 
emission lines, leaving an approximate Fe II emission spec- 
trum. The V04 template is also based on an observed spec- 
trum of 1 Zw 1, but was constructed by removing a much 
larger list of non-Fe II emission features, and also removing 
narrow permitted and forbidden Fe II lines which were found 
to be much stronger in 1 Zw 1 than in most typical Seyfert 1 
spectra. Each Fe II line in the remaining spectrum was then 
modeled with a Lorentzian profile of FWHM= 1 100 km s~', 
producing a noise-free template. The more rigorous removal 
of non-Fe II emission features and narrow Fe II lines by V04 
is a major reason to prefer the V04 template over the BG92 
template for these fits. Furthermore, the empirical and some- 
what subjective procedure used to remove the AGN contin- 
uum by BG92 may be responsible for the fitting degeneracy 
between the Fe II and FC components described above. If the 
AGN continuum were undersubtracted when constructing the 
BG92 template, there would be almost no noticeable impact 
on the quality of individual spectral fits that used the template, 
but the mixing of FC and Fe II emission would likely bias the 
reverberation lags measured for the Fe II component. 

For NGC 4593, the CCF measured between the Fe II and FC 
light curves is certainly better when the V04 template is used, 
compared with the BG92 template. The reason for this differ- 
ence is not clear, and it only appears to occur for NGC 4593 
and not Mrk 1511, but it does add one additional reason to 
prefer the V04 template fits. For all of these reasons, we con- 
sider the measurements done using the V04 template fits to 
be our best-quality results. Most importantly, our main results 
are independent of the template choice: in either case we find 
that Fe II emission does show a strong reverberation response, 
with a lag time longer than that of H/3. As described in the Ap- 
pendix, we find similar results when using the newly released 
multi-component Fe II templates of iKovacevic et al] (120101) . 
With these new templates, we obtain slightly different values 
of Teen for the cmission lines, but the Fe II lags are still consis- 
tently longer than those of H/3, with rcen(Fe II)/rcen(H/3) w 1.5. 
This further confirms that the larger size of the Fe Il-emitting 
region (compared with H/3) is genuine and not an artifact of 
template choice. 

It is also worth noting that the different template fits yield 
differing amounts of Fe II and He I emission underlying the 
broad H/3 profile. The uncertain amount of blending of H/3 
with other emission features (each having a different lag rela- 
tive to the continuum) could represent a limiting factor for the 
accuracy of high-fidelity reverberation measurements, partic- 
ularly for velocity-resolved lag measurements which attempt 
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to determine the distribution of lag across the velocity width 
of the line. Any ambiguity in decomposing the H/3 spectral 
region into different emission-line components should be con- 
sidered as a source of systematic uncertainty in determining 
the H/3 lag distribution across the line profile. 

5. DISCUSSION AND CONCLUSIONS 

This is the first time that such clear reverberation sig- 
nals have been seen for the optical Fe II blends in Seyfert 
galaxies. Our high-cadence monitoring data reveal that the 
Fe II emission in these galaxies does reverberate on short 
timescales in response to continuum variations, with a well- 
defined cross-correlation lag time. This gives direct evidence 
for an origin of the Fe II emission in photoionized gas in the 
BLR. The ambiguit y in previous Fe II reverberation results 
for NGC 5 548 ( Vest ergaard & Peterson 2005) and Ark 120 
(iKuehn et a l. 2008) might be attributable to the lower cadence 
of the monitoring data used in these studies, although in both 
cases the monitoring duration spanned several years. It is 
also possible that the Fe II variability behavior in our two ob- 
jects is not representative of the entire population of Seyferts. 
Measurement of Fe II variability over a broad range of AGN 
properties should be a high priority for future high-cadence 
reverberation-mapping programs. By exploring a broad range 
of luminosities, it might be possible to test whether the optical 
Fe n emission follows a radius-luminosity relationship simi- 
lar to that of the H/3 h ne, with a l uminosity dependence of 
approximately r oc LP'^ (iBentz et al. l 2009). It would be partic- 
ularly interesting to test whether the ratio of Fe II to H/3 radii 
varies systematically along the Eigenvector 1 sequence as a 
function ofL/LEdd- If the Fe Il-emitti ng gas is infalling t oward 
the central engine, as proposed by iHu et al.l (l2008bl) . then 
there might also be observable correlations between the re- 
verberation lag and the redshift of the Fe II lines, although ev- 
idence for radial inflow remains controversial (Sulentic et al. 
l20T2h . 

The measured values of Teen give an approximate mean ra- 
dius for the zone from which each line is emitted within the 
BLR. Comparing the lags of Fe II and H/3 for the preferred 
V04 template fits, we find Tce„(Fe II)/reen(H/3) = 1.9 ±0.6 
and 1.5 ±0.3 in NGC 4593 and Mrk 1511, respectively. This 
gives a direct indication that the Fe II emission arises predomi- 
nantly in the outer portion of the BLR on larger scales than the 
H/3 emission region. This conclusion is consistent with sev- 
eral other recent lines o f evidence for an outer BLR locatio n 
for the Fe II emission ('Sluse et al."200r; 'Hu et al."2008a'b'; 
Matsuoka et al. 2008; PoDovic et al. 2009; Gas kell 2009; 
Kpvacevic et al. 2010; Shields et al. 2010; Shapo valova et al.l 
2012; .Hu et aL.2012;.Mor & Netzer„ 2012). We use the term 



"outer BLR" to denote emission from a region having a larger 
mean size than the H/3 emission zone. Aside from Fe II, the 
only broad emission l ines seen to have la gs longer than that of 
n/are C III] A1909 (Peterson & Wandel 1999) and Ha (e.g., 
iKaspi et aL.2000. ; Bentz et aL.2010 a). Thus, the emission re- 
gions for these lines correspond to the outermost observable 
portion of the BLR. The BG92 template fits give different ra- 
tios of Fe II to H^ lag, but we find that Tcen(Fe II) > Tcen(H/3) 
for either template. 

The longer lag of Fe II relative to H/3 does not mean that 
the H/3 and Fe II emission regions are physically distinct; 
in fact, there must be a very substantial radial overlap be- 
tween them. Recent progress in transfer-function modeling 
for reverberation-mapping data has provided direct illustra- 
tions of the radial extent of the BLR as seen in Balmer lines. 



For Arp 151, the transfer functions show that the H/3 emis- 
sion extends over a broad radial zone with an outer extent 
that is more than twice as large as cTcen, and the Ha-emitting 
region is several times larger in extent than the mean H/3 ra- 
dius (iBentz et al . 2010b). If the ratios of Fe II to H/3 size for 
NGC 4593 and Mrk 1511 are typical, then the Fe Il-emitting 
zone would encompass the outer portion of the H/3-emitting 
zone of the BLR and beyond, probably corresponding roughly 
to the region em itting Ha. This echoes the recent finding 
by iHu et aP (1200 8a) that quasar H/3 profiles generically con- 
tain an intermediate-width component which may be emitted 
from the same region that produces Fe II, while the very broad 
component of H/3 would be emitted from smaller radii. Fur- 
ther support for a link between the Fe II region and the outer 
portion of the Bahner-line emitting zone of the BLR comes 
fro m a new princip al-component analysis of quasar spectra 
by iHu et al.l (l2012h . They demonstrate that one of the pri- 
mary eigenspectra of their quasar sample is essentially the 
sum of two components: the Fe II emission spectrum and 
an intermediate-width core component of the Balmer lines, 
which are readily interpreted as arising from a similar spatial 
region^^ 

Hu eTaD (I2008bh found that Fe II line widths in SDSS 
quasars are typically about 3/4 of the H/3 widths. In order 
to test whether NGC 4593 and Mrk 1511 are consistent with 
this trend, we examine the results of the spectral fits done us- 
ing the BG92 Fe II template, since this closely follows the 
fitting method used by Hu et al. (2008b). We assume that the 
total velocity width of FWHM is given by the quadrature sum 
of the FWHM of the broadening kernel from the fit and the 
FWHM of the line profiles in the template itself which we take 
to be 900 km s"' for consistency with lHuetal] (l2008bl) . We 
also correct the observed widt hs for an instrumen tal broaden- 
ing of FWHM« 315 km s"' dBarth et al.ll20TTal) . Then, we 
obtain Fe II widths of FWHM= 3330 ± 153 km s"' for NGC 
4593 and 3128 ± 143 km s"' for Mrk 1511, where the un- 
certainties are based on the night-to-night scatter in the width 
of the Fe II Gaussian broadening kernel determined by the 
fitting procedure. The broad-H/3 component models for the 
two AGNs have FWHM = 4395 ± 362 km s"' for NGC 4593 
and 4171 ± 137 km s"' for Mrk 1511. In both cases, then, 
FWHM(Fe II)/FWHM(H/3) is almost preci sely 0.75 , closely 
consistent with the average result for the iHu et al.l (l2008bb 
sample. 

However, we also find that the inferred Fe II widths seem 
to be substantially dependent on choice of template. Us- 
ing our results from fitting with the V04 template, we ob- 
tain broader Fe II widths of FWHM = 5044 ± 176 km s"' 
and 4459 ± 1 16 km s"' for NGC 4593 and Mrk 1511, in both 
cases broader than the H/3 FWHM values. This appears to 
be the result of a combination of factors including different 
relative strengths for individual Fe II lines between the two 
templates, and the fact that the V04 template is constructed 
using Lorentzian models fitted to each line. One consequence 
of these Lorentzian profiles is that the FWHM values of the 
template itself and of the Gaussian broadening kernel do not 
simply add in quadrature to give the FWHM of the total line 
profile; the FWHM of the Gaussian-broadened template pro- 
file is broader than the quadrature sum of the template FWHM 
and broadening kernel FWHM. 

Based on the fits done using the BG92 tem- 
plate, we can conclude that these two AGNs have 
FWHM(Fe II)/FWHM(H/3) ratios consistent with the 
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mean value found by iHu et alj (l2008bh for a large SDSS 
sample, so there is no evidence that they are outliers from 
the normal AGN population. However, the absolute deter- 
mination of Fe II width seems to be subject to substantial 
systematic uncertainty due to the details of how the fit is 
performed and choice of template. Further investigation of 
the cause of this discrepancy in Fe II widths for different 
templates is beyond the scope of this paper, but this problem 
is likely to affect all inferences about Fe II profile widths 
measured by template fitting, particularly when the individual 
Fe II features are completely blended into a pseudocontinuum 
as is the case in these two objects. 

In addition to the systematic issues related to choice of Fe II 
template, there are some additional and important caveats to 
note about the measured Fe II lags. The CCF methodology 
only gives a single, simplistic measure of the lag time for 
an emission line, whereas in reality a given spectral line will 
be emitted over a broad range of radii. It might be possi- 
ble to obtain more detailed information on the radial distribu- 
tion of Fe II emissivit y by applying t h e geo metric modeling 
methods described by iPancoast et al.l (12011b or by applying 
techniques to extract the shape of the transfer function from 
the data (Bentz et al. 2010b; Grieret al. 2013). However, the 
Fe II light curves are relatively noisy and higher-quality data 
might be required, perhaps over a longer monitoring dura- 
tion, in order to go beyond the simple cross-correlation de- 
termination of the Fe II lag presented here. Furthermore, our 
spectral fits assume a uniform flux scaling for the entire Fe II 
template, which results in an average lag measurement for 
the entire complex of Fe II blends, but this method is unable 
to explore the possibility of different response times for dif- 
ferent Fe II lines. Recent work has improved on traditional 
template-fitting methods by allowin g for different behavio r 
among different groups of Fe II lines VKovac evic et al.ll2010h . 
and Shapovalova et al. (2012) demonstrated that in Ark 564, 
Fe II lines from different multiplets showed differing levels 
of correlation with continuum variations. Application of such 
methods to reverberation-mapping data could potentially de- 
tect or constrain differences in reverberation timescale for dif- 
ferent Fe II multiplets as well. 

The primary conclusions of this study are that measure- 
ment of the reverberation lag of the optical Fe II blends is in- 
deed possible in favorable cases, and that in these two AGNs 
the Fe II emission responds directly to continuum variations 
with a lag time that corresponds to the outer portion of the 
BLR, somewhat larger than the H/3-emitting radius. With 
well-sampled data it is possible to detect this reverberation 



signature with a level of significance comparable to that of 
high-quality H/3 reverberation measurements. Detection of 
Fe II reverberation does require some special circumstances, 
in particular a high amplitude of continuum variability, as well 
as relatively strong Fe II emission overall. In all other ob- 
jects from our 201 1 sampl e, either the Fe II emission was too 
weak (as in Mrk 50; Bar th et al.ll2011bh . or the overall flux 
variability was too low for measurements like these to be suc- 
cessful. Finally, these results illustrate the value of carrying 
out spectral decompositions of reverberation-mapping data, 
in order to measure accurate light curves for Fe II, He II, and 
other weak or low-amplitude spectral features. Application of 
these methods for measurement of emission-line light curves 
for our entire 201 1 sample will be described in future papers 
in this series. 
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APPENDIX 



After this paper was submitted, the referee informed u s that a public release of the multicomponent Fe II templates from 
iKovacevic et al.l (120101) had recently been announced (Popo vic et al.l l2013i). Although a full study of Fe II variability using these 
new templates is beyond the scope o f this paper, we canied out an initial examination by adapting our fitting method to use these 
templates. As described by lKova5e vic et al. (2010) and Shapovalova et al. (2012), the template set includes spectra describing 
Fe II lines from four separate multiplet groups (denoted as the F, G, S, and P groups), as well as an additional template containing 
other lines found in the 1 Zw 1 spectrum. Carrying out fits using this set of five templates gives substantially more freedom to 
accurately fit observed Fe II spectra than the monolithic templates of BG92 or V04, at the cost of adding four additional free 
parameters to allow for the individual flux scaling of each Fe II component. 

Figure|6]illustrates the fits to the mean spectra using these multicomponent templates. Similar to the fits with the V04 and BG92 
templates, the model reproduces the over all spectral shape well, but there are significant differences in the fit details compared 
with the monolithic templates. Using the IKovacevic et aP (2010) templates, the fit assigns relatively more flux to the starlight 
component and less to the FC. Also, these fits force the fluxes of all three He I lines to zero, and the red shelf of H/3 becomes 
dominated by Fe II, similar to the results from using the BG92 template. Another notable difference is seen at the shortest 
wavelengths, below about Aiest = 4400 A. In this region, the primary Fe II contribution is from the P group, and the fit forces the 
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Fig. 6. — Same as Figure|2] but for fits done using the lKovacevic et aH 120101) Fe II templates. The grey curve represents the sum of the five individual Fe II 
components. 

normalization of this component to zero. This is most likely a spurious result due to degeneracies in the model fitting process, 
but it occurs consistently when fitting each individual spectrum for both AGNs. These multicomponent templat es are optimally 
suited for use with AGN-dominated spectra having very strong Fe II emission, such as the objects studied by Kovacevic et al.l 
( 120101) and lShapovalova et al.l ( 12012) . but the high starlight fraction and relatively weaker Fe II emission in our targets presents a 
more difficult case study for constraining the weights of the five-component Fe II model. 

In order to check for differences with respect to the previous spectral fits, we measured H/3 and Fe II light curves based on these 
decompositions following the same methods described previously, and carried out cross-correlations of these light curves against 
the y-band continuum. The Fe II light curve was integrated over the same spectral region used for the monolithic templates, in 
this case corresponding primarily to flux from the F group template. 

For NGC 4593, we find Tcen(H/3)= 3.92!°;™ days, and Tee„(Fe II)= 5.97![ ;i^' days. For Mrk 1511, the results are Tce„(H/3)= 
5.76;!;q94 days and rcen(Fe II)= 8.90;!;['25 days. These results are generally consistent with the measurements done using the BG92 
and V64 template fits, with the largest disagreements only being at slightly greater than the la level. The modest disagreements 
further highlight the fact that the reverberation results are somewhat sensitive to the different structures of the Fe II templates. 
Despite these differences, our primary result is essentially unchanged: the Fe II reverberation lags based on the new template fits 
are ^ 50% longer than those of H/3, pointing to an origin for the Fe II emission in the outer portion of the BLR. In future work, 
these multicomponent templates may prove to be most advantageous when fitting very high S/N spectra of AGNs having smaller 
starlight contributions and stronger Fe II lines whose relative amplitudes can be more tightly constrained in the fits. 
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